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Abstract. We discuss various compatibility criteria for overde- 
termined systems of PDEs generalizing the approach to formal 
integrability via brackets of differential operators. Then we give 
sufhcient conditions that guarantee that a PDE possessing a Lie 
algebra of symmetries has invariant solutions. Finally we discuss 
models of equations with large symmetry algebras, which eventu- 
ally lead to integration in closed form. 



Introduction 

Overdetermined systems of PDEs always have compatibility condi- 
tions. If these are satisfied, the system is called formally integrable (we 
assume regularity throughout the paper) and we can formally param- 
etrize the space of solutions [C3 IBCG^I [KL^ . 

For Frobenius type systems the compatibility conditions are just 
equalities of mixed derivatives. For instance (we denote the partial 
derivatives as usual by indices) the system £ on M^(a;^,x^) 

{Fij : Uij = fij{x, u, du) \ = 1, 2} 

has compatibility conditions 



Vi{F22) = V^iFu) ^ 1 I?i(/22) = V^if: 



mods. 



112) 



In general this is wrong. For example, consider the following system S 



on M^fx^, x^, x^, x^) 




Ul3 = U22 + fl{x,U,du), 

Ui2 = Uo3 + f2{x,u,du), (1) 
Uo2 = un + f3{x,u,du). 

Equalities of the mixed derivatives do not yield compatibility here, 
because the first compatibility condition involving only two equations 
is of order 3 (the corresponding syzygy operator, which is the Mayer 
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bracket described in Theorem [21 is of the second order), and it depends 
on derivatives of / up to the 2nd order (for more details see Section [T]). 

There exists however a compatibihty condition of order 2 (the syzygy 
operator involves only first derivatives of /; Fi below denotes the dif- 
ference of the left and right hand sides in the respective equation): 

I?i(Fi) + V2{F2) + V^{F^) = mod S. (2) 

In Section [T] we discuss how to calculate the compatibility conditions. 
This is related in the next Section to classical and higher symmetries, 
and we discuss the question of existence of invariant solutions. 

Consider a PDE £ and a subalgebra Q C sjm{£); at this point we 
restrict to the local extrinsic (point, contact or higher) symmetries. In 
general it is not true that S has ^-invariant solutions, as was noticed 
in |Ki| . The simplest counter-example constitute linear systems with 
the symmetry being shift by a nonzero solution. 

Another example is given by a pair of constant coefficients non- 
homogeneous linear scalar PDEs. They commute provided there are 
no zero order terms, and so the equations are symmetries of each other. 
Generically the PDEs are compatible, but this does not happen always 
(e.g. the system {u^x = ^,Ux = 0} has no solutions at all). A similar 
story happens for matrix differential operators and symmetries. 

In Section [2] we demonstrate that this is a consequence of either de- 
generacy or higher dimensionality of the symmetry algebra compared to 
the amount of independent variables. We will prove that under certain 
genericity conditions the symmetry is compatible with the equation. 
This yields (a similar result is proven by another approach in [I V] ) : 

Provided the symbols of S and Q are generic, the system 
£ has Q -invariant solutions. 

In Section [3] we discuss implications that existence of a large sym- 
metry group has on the solution space of a PDE system. At this point 
we need to consider intrinsic symmetries. Relations between extrinsic 
and intrinsic symmetries are given by Lie-Backlund type theorems, see 
[AH IKLV[ IAKO[ lAKj , and this relates this integrability problem with 
what is discussed above (see also [L7], [K2i IGal IBCAj for applications of 
symmetries and generalizations). 

We will concentrate on exact solvability of ODEs and PDEs and 
discuss relations with Darboux integrability via examples. Informally 
a lot of symmetry implies exact integrability, more precisely this holds 
true for maximal symmetric models (can fail for sub- maximal cases). 

We will briefly discuss some symmetric models. Relations of inte- 
grability to transformations and differential substitutions in PDEs is 
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central in IL2I IK3I IK4I for the case when the system £ depends on 1 
function of 1 argument (so called Lie class u = 1). This applies in the 
other cases too, but will be considered elsewhere. 

The paper is organized so that all sections can be read independently. 

Acknowledgement. I thank V. Lychagin and A. Prasolov for 
helpful discussions. Some calculations in Section [3] were performed 
with Maple package DifferentialGeometry by I. Anderson. I 
am grateful for organizers of the conference SPT-2011 for hospitality. 

1. Compatibility, differential syzygies and brackets 

An overdetermined system of differential equations can be viewed 
geometrically as a finite sequence of submanifolds £k C J'^(vr) in jets 
with E^^li D Sk) where vr : — )■ M is a vector bundle and E^^li denotes 
the prolongation (locus of the derivatives of functions specifying £k~i)- 

Formal integrability is equivalent to the claim that all projections 

TTfc+i k '■ £k+i — ^ £k are submersions. Then we can define £^00 = \vai£k- 

' <— 

Let T be the model tangent space for independent variables and 
the model tangent space for dependent variables. Define the symbols 
of order k as Qk = 'i^ei{dTXk,k~i '■ T£k ~^ T£k-i) C S'^T* ® and 
let g = ®gk be the symbol bundle over £. For k greater than the 
(maximal) order / of £^ we can define gk = S^'^T* ® giH S''T* (g) N. 

The Spencer (5-cohomology group H*'*{£) is the cohomology of the 
complex (yf* ® A*(T*) with symbolic de Rham differential over £ [Sj. 

The formal theory of differential equations identifies obstructions 
to formal integrability (compatibility conditions) as certain elements 
Wk G H^~^''^{£) (structural functions or Weyl tensors). For geometric 
problems H*''^{£) is the space of curvature/torsion tensors IKL3 . 

Another way to look at this is to regard compatibility conditions as 
differential syzygies, i.e. relations between generators of £. 

Denote generators of this system by F = {Fj}^^^, so that £ = 
{Fi[u] = 0, . . . , Fr[u] = 0}. Then the symbol spaces can be expressed 
through the symbols of linearizations of these (nonlinear) differential 
operators (we will use the notation a for the symbol) 

g = KeT{a{iF,), a(£f J} C ST* ® N. 

It turns out that g* = ®gl is an ST-module, where the latter is 
viewed as the algebra of polynomials on T*. This is called the symbolic 
module of the system g* = M.e (also known as the symbol module [SJ). 

Definition 1. A differential operator G from the left differential ideal 
(F) is called a differential syzygy for the system £ if its symbol (j{(!.g) 
is an ordinary (algebraic) syzygy for the symbolic module M.s- 
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If £ is linear we can consider (F) as the left module over the algebra 
Diff(l, 1) of scalar linear differential operators (on the trivial rank 1 
bundle 1 over M). In the nonlinear case, this should be changed to 
the algebra ^ Diff(l , 1) = C°°(J°°7r)®c-(M)Difr(l, 1) of ^-differential 
operators, see jKl^ : (F) = (A o | A e ^Diff(l, 1)). 

Thus to every differential syzygy there corresponds an algebraic 
syzygy, but because the operator of taking symbol is not injective, 
the reverse problem 

syzygy -w differential syzygy, 

is not uniquely solvable. Moreover it is not possible to construct q as a 
homomorphism, but we would like to have a map q, defined on a finite 
generating set of syzygies, with nice algebraic properties (extension to 
the whole syzygy module is similar to quantization theories). This will 
be done below for the class of (generalized) complete intersections. 

Let fi = cr(£i?.) G ST ® A^* be symbols of the differential operators 
defining £ and let J^difi = be a syzygy, with gi G ST being some 
polynomials on T* with ki = deg(/j) = ord(Fj). 

Choose any "^-differential operators Gi with a{Gi) = gi. Then the 
corresponding differential syzygy has the form 

V = ^G,oF,. 

The order of this operator is — 1, where k = maXj{ord(G'j) -|- ki} is 
called the order of the syzygy. 

Non- uniqueness comes through the lower order terms in Gj. So the 
class V modj7fc_i(F) is well defined, where 

Jt{F) = {J2Qi°Fi- ord(Qi) < t - h}. 

Denote by [S] the equivalence class S mod J7a:-i(F) of the differential 
syzygy S = q{s), where k is order of the algebraic syzygy s (notice that 
this equivalence class is independent of the choice of q so far it is right- 
inverse to the symbol map a). 

Theorem 1 ( IKL4I ). An overdetermined PDE system £ is formally in- 
tegrable iff for a basis {sj} of algebraic syzygy the corresponding classes 
of differential syzygies vanish [Sj] = 0. 

This follows from the fact that the Spencer (5-cohomology group 
ffk-i,2^g^ equals the corresponding graded second cohomology group 
for projective resolution of the symbolic module Ms, i.e. it can be 
enumerated via a basis of algebraic syzygies. 

There are two basic approaches to construct the arrow q: 



SYMMETRY, COMPATIBILITY AND EXACT SOLUTIONS 5 

1. Construct differential syzygies successively in order k, i.e. ac- 
cording to passage Sk Sk+i- This corresponds to prolongation- 
projection approach having origin in E. Cartan's equivalence method. 

2. Successive identification of differential syzygies involving two op- 
erators Fi,Fj, then involving three etc. This represents compatibility 
as Massey products IKL4I and is related to deformation of the symbolic 
module Ai^ (or to the corresponding noncommutative D- module S*). 

We will elaborate the first idea for PDE systems with nice character- 
istic variety ChaT^{£) = G P^T* : rank[(T(£i7)(0] < m = dimA^}. 

Now I want to present the explicit form of compatibility conditions 
for certain generic overdetermined systems of PDEs. We start with the 
scalar case m = 1 (single u); n = dimM is arbitrary. 

The Jacobi bracket {F, G} of scalar (non-linear) differential opera- 
tors F,G & diff (1, 1) is defined via the linearization operator as follows: 

{F,G}=iFG-iGF 

If ord F = k,oTdG = I, then ord{F, G} = k + l-l. 

Let {x\Ua), 1 < i < n, a = (ii, . . . be the "canonical" coordi- 
nates on the jet-space, Vi = d^i + Y^ Ua+udu^ the total derivative w.r.t. 
X* and let = . . . P^" be the operator of higher total derivative. 
We can express the Jacobi bracket as 

= Y.V^{F)duAG) - V„{G)duAF)- 

We define the Mayer bracket of Fj, Fj by the formula 

[Fi,Fj\e = {Fi,Fj} mod Jfc,+fc^,_i(F). 

For first order scalar operators these are respectively the classical La- 
grange and Mayer brackets. 

Theorem 2 ( |KLi| ). Consider a scalar system S C J^{M) given by 
r <n differential equations Fi[u] = 0, . . . , Fr[u\ = 0, such that for each 
point Xk & £ the characteristic varieties for the equations Fi are jointly 
transversal, i.e. codim[Char!^^(£) <z¥ {T*M)'^^^ =r. 

Then the system is formally integrable iff all the Mayer brackets van- 
ish: 

[F,,F^]e = 0, l<i<3<r. 

The Koszul complex is the minimal resolution of the symbolic mod- 
ule Ais for complete intersections, whence the algebraic syzygies are 
generated by commutators. Thus the arrow q associates the higher 
Jacobi bracket {, } to the commutator [, ] (and so can be treated as a 
quantization) . 
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If the condition of complete intersection is violated, then the con- 
clusion of the theorem could be wrong. To see this let us re-visit the 
example from the Introduction. 

Example 1. For system ([1]) the compatibility conditions [Fi, F2]£ = 0, 
[F2, F-sls = 0, [F3, Fi]^: = should necessarily hold (Fj denotes the 
difference of the left and right hand sides in the i-th equation of ([T])), 
but they do not form a basis in differential syzygy module (and so are 
not sufficient for compatibility). A basis is given by relation and a 
similar relation 

Vo{Fi) + I?i(F2) + V2{F3) = mod £. 

Bracket-relations are differential corollaries of these two, which yield 
a basis in the second cohomology group (counting compatibility con- 
ditions) H^''^{S) = R^. It is easy to see that the system S is not a 
complete intersection: the characteristic variety is the normal cubic 

Char^(£:) = G CP^ I ^1^3 = ^2, = ^06, ^06 = 
= {[A^: A^iA: 1]|AgC = CUcx)} 

and it cannot be given by 2 = codim Char^(£^) equations. 

We can however have other explicit formulae for scalar non-complete 
intersections. Let us point out one example. The following statement 
can be proved similarly to Theorem [2] (and it follows from Theorem [1]). 

Theorem 3. Consider the system 

S = {FioGj + ... =0\i = l...r,j = l...s} 

where dots stand for the lower order terms. It has reducible character- 
istic variety Char'^(£) = {a{£pj = 0} U {o"(£g'J = 0}, and we suppose 
that the intersection {(r{ip.) = 0}n{a{iGj) = 0} has codimension r + s . 
Then if ki = ord(-Fj), Ij = OTd{Gj), the compatibility conditions are 

{Fi,F,}Gk = modJk^+k,+i,-i{£),Fi{Gj,Gk} = modJk,+i^+i,-iiS). 

Let us now consider vector systems of PDEs on m = {u^, . . . , u"^). 

Definition 2. A system S C J''{n) of PDEs = 0, . . . , Fju] = 
is called a generalized complete intersection if 

(1) m<r<n-\-m — 1, where n = dimM, m = rank(7r); 

(2) The complex projective characteristic variety Char^(£^) C P^T* 
has codimension r — m + 1; 

(3) The characteristic sheaf K, over Char'^(£') has fibers of dimen- 
sion 1 everywhere. 
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Let US introduce a multi-bracket of linear (scalar) differential opera- 
tors Vj G Diff(m • 1, 1) by the formula 



rra+l 

7jr"| l<i<m 

fc=l 



{Vi, . . . , = Ndet ■ 



where Ndet is a version of non-commutative determinant. For non- 
linear differential operators Fi G diff (vr, 1) the multi-bracket is 

{Fi, . . ., Fm+i} = 

agSm,/3GSm+i 

Then we define the reduced bracket 

Theorem 4 ( IKL4I ). A system of generalized complete intersection 
type 

S c J^ir) = {F^[u] = 0, . . . , Fr[u] = 0} 
is formally integrable iff all the multi- brackets vanish due to the system: 

[Fjj, . . .,Fi^^^]£ = 0. 

For such systems S the map q associates the multi-brackets of dif- 
ferential operators to certain determinental identities, and these multi- 
brackets satisfy the Pliicker identity |KL^ 



2. Symmetries and compatibility 

Consider a compatible system £ = {Fi = 0, . . . , = 0}. Let Q be 
a subalgebra of the algebra sym(£) of classical or higher symmetries 
|KLV] . written as another PDF system {Si = 0, . . . , 5"^ = 0}. 

Consider at first the scalar case, when m is a single function. Denote 
ki = ord(-Fj), Ij = ord(S'j). Then the symmetry condition is 

{F„5J = 0modJfe^+i^_i(^). (3) 

In addition we have from the symmetry condition that 

{S,,S,} = mod Ji^^i^_i{g). (4) 

Theorem 5. Assume that the joint system £ + Q: {Fi = 0, Sj = 0} is 
a complete intersection, i.e. its characteristic variety has codimension 
r + k < n. Then this system is compatible. 



For k = 1 this result coincides with Theorem 17 from KLi . 
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Proof. Since the system is a complete intersection, Theorem [2] imphes 
compatibihty provided that the Mayer brackets (on the joint system 
S + G) vanish. Compatibihty of £ yields 

{Fi,Fj} =0 mod Jk^+k,-i{S). 

Vanishing of the other brackets is given by ([3]) and (HI). □ 

The condition of complete intersection is important. 

Example 2. For the KdV equation ut = uUx+Uxxx the algebra of clas- 
sical symmetries is generated by: Tq = ut, Ti = Ux, R = 3tut+xUx+2u, 
r = tUx + 1. For every symmetry (one-dimensional subalgebra) there 
are invariant solutions. Moreover the joint system KdV+ symmetry is 
compatible, so the space of invariant solutions had the expected dimen- 
sion (1 for a symmetry S G (Ti, F) and 3 for any other symmetry). 

However there are no nontrivial invariant solutions for the following 
two-dimensional subalgebras of symmetries: 1) {Tq,R); 2) {Ti,R); 3) 
(Ti, F); 4) (F, R) (only zero solution for the joint system KdV-t- 2 sym- 
metries in the cases 1 & 2 and nothing at all in 3 & 4). And indeed, 
in this case the joint system £^ -|- ^ is not a complete intersection. 

A similar observation on (non-)existence of invariant solutions ap- 
plies to the infinite-dimensional higher symmetry algebra of the KdV. 

Now let us discuss the general case, when the unknown -u is a vector- 
function. Consider at first the following example of matrix linear dif- 
ferential operators: 

/ VxV, 0\ (VxVy + Vx 

[VxVy-Vx Vx) V ^-^1/ 

They commute [A, 5] = and so do the inhomogeneous operators 

^-^•(:)-(o)- - 

In other words G G sym(F). However the operators are incompatible. 

We shall show that for two generic nonlinear differential operators 
the condition G G sym(F) implies compatibility of the joint system, 
and so existence of invariant solutions (in the amount given by the 
usual formal calculus of dimensions IC3I ). 

This result was noticed by S.Igonin and A.Verbovetsky and a proof 
using a different method will appear in [IV]. We would like to give an 
explicit criterion for this compatibility. At first we notice the following 

Lemma . Consider a system £ of r > m differential equations for m 
unknown functions. Its characteristic variety Char''^(£^) has codimen- 
sion < r ~ m + 1. 
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Notice that co dimension cannot exceed n since we adopt the conven- 
tion dim0 = — 1 and the ambient space satisfies dimP'^T* = n — 1. 

Proof. The symbol of the system is a r x m matrix P with polynomial 
entries, and the characteristic variety is 



Denote by Ajj...j^ the determinant of the m x m minor generated by 
rows ii, . . . , im- Char'^(£^) is given by the conditions Aj^...j^(^) = for 
alHi < ■ • • < im, but this collection is excessive. 

We can suppose that Ai ^ ^ (if all the minors are degenerate 
there is nothing to prove). For j > m let us denote by Pj the upper- 
row submatrix of P of size j x m. By induction we can suppose that the 
set Sj = G P^T* : rank(Pj (^)) < m} has codimension < j — m + 1. 
Consider the subset where the above rank is < m — 1 . 

If Tj'j has a component in S^, then addition of a row to Pj does not 
increase the rank over m — 1, and the codimension of S^+i containing 
this component is the same as for Sj. 

Otherwise S^- C Sj has codimension at least 1 and in the complement 
the matrix Pj{C,) has rank (m — 1). Near every point ^ we can choose 
(m — 1) X m subminor of maximal rank. Adding to it the row number 
(j -|- 1) we get a m X m matrix whose determinant we write as Aj+i. 
Then the defining relation for Sj+i C outside is A^+i = 0, 
whence the relative codimension is 1. 

Thus codimension of S^+i in the complex projective variety is at 
most j — m and this gives the induction step. □ 

In particular, for r = 2m the codimension is at most (and generically 
for m < n it is) m + 1. We shall refine this in the case, when the system 
£^ is a PDE coupled with its symmetry. 

Consider the PDE system F = G = 0, where F, G G diff (tt, tt) 
are nonlinear differential operators from a rank n bundle tt to itself 
and G G sym(F). Denote the symbols of these operators by P = 
(j(£i?), Q = o"(£g)- The symmetry condition {F,G} = modF implies 
the following relation on polynomial matrices P, Q with some other 
polynomial matrix K 



Proposition 6. For m > 1 let P and Q be two homogeneous polyno- 
mial m X m matrices satisfying ([^. Then the characteristic variety 



1^ g pL^-* . rank(P(0) < m}. 



PQ = KP. 



(5) 
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(characteristic variety of S : {F = G = 0}) has codimension < m. 

Proof. Denoting by Ri{A) the row of matrix A, we have from the sym- 
metry condition: 

m 

J2mRk{Q) ^ {Ri{P),...,RUP)) V« = i,...,m. (6) 

k=l 

We can suppose that detP(^) ^ (otherwise the claim follows from 
the Lemma), so that the equation det P{C,) = determines a subvariety 
S C W^T* of codimension 1. 

Let So C S be given by the equation P{^) = 0. First let us study 
the points C, E T,\ Eq. At such ^ there exists an entry ^ 0. Then we 
write the conditions Rj{Q) G (i?i(P), . . . , Rm{P)) for j ^ k. These are 
no more than m — 1 equations and so they specify a subvariety /C C S 
in a neighborhood of ^ on which also RkiQ) G (i?i(P), . . . , Rm{P)) by 
([6j). Thus this /C is a part of Char''' of codimension < m. 

Now let us consider Sq. If it has codimension > m, it is negligible 
or is a part of Char*^ by the above argument. But if its codimension 
is < m, then we must add the condition det (5(0 = specifying a 
subvariety /Co C So of codimension < 1. Since P(0 = for ^ G So the 



rank of the matrix 



P{0 

codimension < m. □ 



is < m. Thus this /Co is a part of Char^ of 



This Proposition is important, so we would like to indicate an idea 
behind an alternative proof. It will be shown later that condition ([2D 
can be changed to [P, Q] = without loss of generality, so we adopt 
this condition. 

By Gerstenhaber theorem |Gej every pair of numeric (complex) com- 
muting matrices P, Q is contained in a commutative algebra (with 1) 
of dimension m. If one of them has simple spectrum, say P, then by 
Cayley-Hamilton theorem this algebra is generated by {P*}™q^ fZ]. 

This being generalized to matrices with polynomial entries, would 
imply existence of polynomial matrices Zq, Zi, . . . , Z^-i and the scalar 
polynomials aj, bj such that 

P = aoZo + aiZi + - ■ ■ + am-iZm~i, Q = boZo + biZi + ■ ■ ■ + bm-iZm~i- 

Therefore one of the strata of the characteristic variety is given by 
the condition that the vector (ao, . . . , am-i) is parallel to {bo, . . . , 6m-i) 
(m — 1 equations) plus one equation det P = (or det (5 = 0), implying 
codim Char*^ < m. 

Let us demonstrate this in the case m = 2. Since identity matrices 
commute with everything, we can subtract a multiple of them to make 
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'fPu 


Pl2\ 


( Qu 


qi2\ 


-\P21 


J 


U21 


\ 



P22 = <i22 = 0. Then commutativity [P, Q] = yields: 

^ P}]_ ^ P}1 ^ P2i 
qu qi2 q2i ' 

which means that these matrices are proportional to some polynomial 
matrix. Indeed, the first equality gives pu = riSi, qu = riS2, Pu = 
^^251, qi2 = ^282 for some polynomials r^, Sj. 

Substituting this into the second equality we get refining: si = uiVi, 
S2 = U1V2, P21 = U2V1, q2i = U2V2, whence the claim P = agl + aiZ, 
Q = bol + hZ. 

Let us return to sufficient conditions for compatibility of symmetries. 

From what was shown above it follows that in the non-scalar case 
S = {F = G = 0} is never a generalized complete intersection (so 
compatibility cannot be deduced on the basis of Theorem H]). This 
is a consequence of non- commutativity of the matrix algebra; another 
feature of the algebra of matrix differential operators (m > 1) is that 
generically ord{F, G} = ord(-F) + ord(G). 

Theorem 7. Let F = be a determined PDF and G = its symmetry 
(both mxm systems). If for the joint system S: F = G = the variety 
Char''(£^) has codimension m, then £ is compatible. 

The idea is that since this is the maximal possible codimension, there 
is no space for another syzygy except for the symmetry relation. 

Proof. Let us first give the proof in a particular case, when both F and 
G are of the first order. Since the equation F = is determined (det P 
is not identically zero), we can write it in the evolutionary form. 

Thus F[u] = ut — Fo[u], where Fq does not involve Ut terms but 
contains some other derivatives of m (m is a vector function with m 
components). We can substitute Ut = Fq into G = and get a PDE 
that is free of Ut terms. We continue to write G for this new operator. 

Denoting the symbols of F and G by P and Q, the symmetry condi- 
tion is PQ = KP for some polynomial matrix K. Since P = t1 + Pq, 
where r stands for the symbol of Vt differentiation and Pq (as well as 
Q) is free of r, we conclude that K = Q, i.e. [P, Q] = 0. 

We claim that G is also determined, i.e. det Q ^ 0. Otherwise the 
characteristic equation is given by m — 1 equations at generic point. 
Indeed, the proof of Lemma/Proposition [6] can be rewritten to start 
with equation det Q = 0, which is identically true and so the number 

of defining relations for the variety : rank < m} is by 1 less 

than in the general case. This yields codim Char'^(£^) < m and so 
contradicts the assumptions. 
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Thus the system G = is compatible (as well as the system F = 0). 
Due to evolutionary form of P (and absence of r in Q) there are no 
syzygies of the system P = Q = besides the commutation relation. 
So the claim follows from Theorem [H 

Now let us consider the general case, when the orders can exceed 
1. Let us re- write the system in the 1st order evolutionary form. We 
keep the same letters F and G for the differential operators. This 
invertible transformation does not affect the symmetry or compatibility 
properties and it preserves Char^(£^), but it increases m. 

This increase is composed of introduction of new variables and re- 
writing the operators F and G via them. The first type equations 
(including equalities of mixed derivatives) have no influence on the 
characteristic variety and when coupled with new F they give the same 
codimension 1 sub-variety in CP""^ as the old F. The new equation G 
(re-written through the new variables and not involving t-derivatives) 
still consists of m equations, and by the same argument as above it is de- 
termined. The symmetry condition written as 'Dt{G)\ut=Fo[u] = (now 
u consists of old and new dependent functions) is the only differential 
syzygy, and we conclude that S = {F = G = 0} is compatible. □ 

Example 3. Consider the Kadomtsev-Pogutse equation (which is a 
simplification of the MHD system): 

F2 = ^xxt + 'fyyt + 'fx'fxxy + ^x^yyy " ^y'f XXX ^y'-fxyy 

'^xxz '4^yyz '^x'^xxy '^x'^yyy ~l~ '^y'^xxx ~l~ '^y'^xyy 0. 

Its symmetries consist of 7 functional families and 2 additional fields, 
see |KV] . For instance here is the first family: 

Ai{a) = a' ■{(f + ^) + a-{Lpz + ^Pt), 

A(a) = a' ■ {ip + tp) + a ■ {iJz + A), 

where a = a{z + t). The characteristic variety of the system £ = {F = 
^ = 0} has codimension m = 2, 

Char^(^) = {[p, -.py-.pz-. pt] eCP'' : pz+pt = 0k 

[Px + = V {pz + IpxPy - IpyPxf - {Pt + VxPy - VyPxf = ]}. 

Thus by Theorem [7| the system £ is compatible, and so gives rise to in- 
variant solutions (see [KVj and references therein for the exact solutions 
of the Kadomtsev-Pogutse equation arising via symmetry methods). 

As we see from the above proof it is not necessary to require that 
codim Char'^(£) = m. What we want to achieve is that the operator G, 
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when F is written in evolutionary form and Vt derivatives are removed 
from G, is determined. This is clearly a generic property. 

Let us now give a sufficient condition in a special case when the 
symbol P = (j{F) (m x m matrix with polynomial entries) can be 
diagonalized via an invertible transformation over polynomials. As 
will be seen from the proof, codim Char'^(£) = 2 and so this case does 
not follow from the above theorem (if m > 2). 

In the proof we will calculate the Matmxm('S'T)-module 

Syz(P, Q) = {{A, B) G Mat^x2™(5T) : AP + BQ = 0}. 

It encodes all syzygies which are the symbols of differential syzygies, 
and the latter are the compatibility conditions by Theorem [TJ 

Theorem 8. Assume that P = cr{iF) is diagonalizable over the al- 
gebra ST. Assume also that the characteristic variety of F = has 
codimension 1 (i.e. F is determined) and has no multiple components. 

Let G G sym(F) he a symmetry (both F, G are mxm PDF systems). 
If for the joint system £ : F = G = the characteristic variety has 
codimension > 1, then S is compatible. 

Proof. We can assume from the beginning that P is diagonal P = 
diag(pi, . . . for some polynomials on T*: pi G ST. Denote the 
symbol of the symmetry G by Q. This is also a matrix with polynomial 
entries, Q = [Qij]inxm- 

The symmetry condition gives the following matrix syzygy: 

PQ = KP QijPi = kijPj (no summation). 

Because the characteristic variety has no multiple components, we con- 
clude from here qij = VijPj, kij = VijPi for i ^ j. This means that 
Q = RP + D and K = PR + D for some diagonal matrix D. Since 
codim Char'^(£^) > 2, pi and qa have no common factors. 

Then the syzygy AP + BQ = is {A + BR)P + BD = 0, which yields 
B = -LP and A = LD-BR = LK. Consequently the MaXmxm{,ST)- 
module Syz(P, Q) is 1-dimensional and generated by {K,—P). The 
result follows from Theorem [H □ 

Remark 1. Based on the generators of the syzygy module Syz(P, Q) 
one can check that multi-brackets of the rows of the 2m x m matrix 

differential operator ^ vanish as a consequence of the symmetry con- 
dition. 



A similar argument shows that an analog of Theorem |5] holds for vec- 
tor nonlinear differential operators, under some genericity assumptions 
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on the characteristic variety/ideal (implying that dimension n shall be 
sufficiently big compared to the number of equations in £ and m > 1). 



n"i'-;,-E^(H-|S)-''. 



Example 4. Consider the non-linear wave equation 

where f * G C°°{M}^"'{t, x^, . . . , x")), i = 1, . . . ,m. Its symbol is diago- 
nal, and the scaling symmetry G[v]^ = C,{v^) for 

n m 

^ = tdt + {l~ 2a) J2 ^^^x^ + 2 XI ^'^''^ 
i=i i=i 

satisfies the property that the characteristic variety of the joint system 
£ : F = G = Q has codimension 2. This system £ is compatible and 
there are invariant solutions. 

3. Exact solvability of differential equations 

According to Cartan's test solutions to a compatible system of differ- 
ential equations £ are formally parametrized by p functions of g argu- 
ments (and some functions with fewer arguments). These numbers are 
expressed through the characteristic variety and the characteristic sheaf 
/C o ver it a s follows: g = dime Char^(£^) + 1, p = deg Char^(£^) ■ dim /C, 
see [BCG3,[KL2]. Denote the abstract space of such functions by 



Closed form solutions refer to parametrization of the generic stratum 
of the solutions space by a differential operator S* : ^ Sol(£^). By 
IC2I for underdetermined ODEs {g = 1) this is tantamounlQ to inter- 
nal equivalence of the equation equipped with the contact distribution 
{£,€£) to some mixed jet space J°"(M, M^). Here a = (ii, ... ,ip) is a 
multi-index characterizing the orders of the dependent variables. 

Example 5. Consider the equation for null curves in Lorentzian space 
of 1 -|- 2 dimensions 

x'itf + y\tf = 1 ^ dt^ - dx"^ - dy^ = 0. 

An obvious solution involves one arbitrary function and the quadrature: 
X = f cos (f)(t)dt, y = J sin (f)(t)dt. But it can also be integrated in the 
closed form: 

t = o-"{t) — cr(r), X = cr"(r) COST + cr'(r) sinr, 

y = —a"{r) sinr -|- cr'(r) cost. 



^In loc. cited only the classical Monge case p = \ [1 equation on 2 unknowns) 
was treated, but the general case is similar. 
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This equation in 1 + n dimensions cannot be solved in closed form 
for n > 2. Indeed internally 

n 
1=1 

The Cartan distribution on £ is generated by Vt = dt + ^idxi and 
the subspace U = TS""-^ = {d^>^ , . . . , ^ x[d^:, =0): Ce = M^V H. 

Let Ai = Cs and A2 = [Ai, Ai] be its derived distribution. Then 
dim Ai = n, dim A2 = 2n — 1 and we have A2 = Ai + {x'^dx^ — x'jdxi)- 
The next derived distribution is A3 = [Ai, A2] = T£. 

Thus by dimensional reasons the only possible corresponding jet- 
space is J'^(R, R"~^), 0" = (1, . . . , 1, 2), with coordinates t,Ui, . . . , m^-i, 



u[, . . . , u'^_i, Its Cartan distribution is 



Ai = {Vt, a^, . . . , where Vt = dt+^u[d^^+u'n-id^'^_^. 

i=l 

The derived distribution is equal to A2 = Ai + . . . , du„_2, 
and A3 = TJ''. 

Though the dimensions coincide, the distributions on £ and J'^ are 
not equivalent. The reason is that 11, which is the maximal involutive 
space of the bracket A^Ai — t- A2/A1, is not the Cauchy characteristic 
space for A2. But in the second case II = {0^'^, . . . , du'^_^, which 
is the maximal involutive space of the bracket A^Ai — t- A2/A1, is the 
Cauchy characteristic space for A2. 

Example 6. Consider the Monge equation w'{x) = {z'{x))'^. The 
general solution depends on 1 function of 1 variable and the form via 
quadrature is obvious, but here is the closed form solution: 

X = a"(r), w = T^a"{T) - 2rcT'(r) + 2a(r), z = rcT"(r) - cT'(r). 

The explanation behind this is the Engel normal form for rank 2 dis- 
tributions in M^. However the next candidate - the Hilbert-Cartan 
equation 

w' = (z'r (7) 

is no longer integrable in closed form (without quadratures) as was 
demonstrated in 1912 by D.Hilbert. In 1914 E. Cartan gave a criterion 



for resolution of underdetermined ODEs in closed form IC2I , which we 
referred to above. 

Already in 1910 Cartan found that the symmetry group of ([7]) is G2 



(though it was not written like this in Ci , he surely knew this) and 
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proved that it is the most symmetric equation among such Monge equa- 
tions with finite-dimensional symmetry groups (linearizable equations 
have infinite-dimensional group of symmetries). 

The more general problem when a PDE is integrable in closed form 
via solutions of a simpler equation (usually ODE) is known as the 
method of Darboux (that is we allow for the operator 5* above to involve 
quadratures and some other nonlocalities). For instance, the Liouville 
equation is Darboux integrable 

„ 2f'{x)g'{y) 
u^y = e =^ u = log . . . , y ^ , 

[fix) +g{y)y 

while the sin-Gordon equation u^y = sinw is not. 

Example 7. The following overdetermined system of PDE on M? ap- 



peared in |Ci| (A is a parameter to be excluded): 

^xx 3A , Uxy — ' ^yy — (^) 

It is a compatible involutive system. The general solution is parametrized 
by 1 function of 1 argument 

X = X, y = z"{t) + XT, 

u = xz{t) + z\t)z"{t) - \w{t) - \rz"{Tf - \t^xz"{t) - |rV, 

where w'{t) = [z'^t))"^. In fact, this system has a common character- 
istic dx — Xdy, which lifts to the Cauchy characteristic of the Cartan 
distribution (of rank 3). The quotient by the Cauchy characteristic is 
a 5-dimensional manifold with rank 2 distribution equivalent to ([7]). 

By a Lie-Backlund type theorem |Ci| the contact symmetry group 
of ([8]) coincides with the internal symmetry group of ([7]), and so is G2. 

According to Goursat [Gou] the general form of overdetermined invo- 
lutive (in this case: compatible with a common characteristic) system 
of 2nd order PDE on the plane is 

r + 2As + A^t = 2^, s + Xt = iPx, t = V'aa, (9) 

where we use the classical notations r = Uxx, s = Uxy,t = Uyy and 
suppose 7^ (nonlinearity). System ([8]) correspondqj to il) = A^/3! 

Removing the last equation from this system we obtain a determined 
parabolic PDE £ of the 2nd order. It has the largest contact symme- 
try group (among non- linear equations) for if) = A^/3! in which case 
excluding A we obtain the Goursat equation 

4:{2s - + {3r - 6st + 2t^f = 0. (10) 



One also need to change y ^ —y to match the sign. 
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This equation has the same symmetry group G2 and it can be parame- 
trized as the 2D tangent cone p(A) + p'{X) to the twisted cubic p(A): 

r = + A^/i, s = iA^ + A/x, t = X + fx. 

Excluding A and /i we get equation flTU]) . 

It has the following geometric reduction to (171), giving exact solutions 
for (fTOl) . The double characteristic dx — Xdy lifts to 

e = P.. - XVy = dx - Xdy +{P- Xq)du - |A=^9p - \X^d,, 

which together with = X^dj. + Xds + dt forms the integrable charac- 
teristic rank 2 distribution Tl (Z Cs. Quotient by it maps the Cartan 
distribution Cs (of rank 4) to a rank 3 distribution A on a 5-dimensional 
manifold M^. This rank 3 distribution is the derivative distribution of 
the unique rank 2 distribution A which maps to zero under restriction 
of the natural bracket A^A TM/A. 

If we identify with M^(?/, u, p, q, A), then A = {dy + qdu + \X^dp + 
Xdq, dx) and A = A + {dg+Xdp). Thus we see that (M^, A) is equivalent 
to Hilbert-Cartan equation ([7j). 

The following deformation of the Goursat equation was studied in 

(4 + e)(2s - ty + (3r - 6st + 2t^f = 0. (11) 

Here e > and for every such number (fTTl) is hyperbolic!! and it has 
maximal symmetry algebra of dimension 9 among all hyperbolic 2nd 
order PDEs on the plane, which are neither of Monge- Ampere nor of 
Goursat type (see [T] for details). The Lie algebra of symmetries has 
Levi decomposition g = SI2 xt (radical r is 6-dimensional) . 

Except for this family there is one more hyperbolic equation with 
9-dimensional symmetry algebra {q has the same form but different r) 

3rt^ + 1 = 0. (12) 

All these equations are Darboux integrable. The last one was studied 
in Goursat |Gou] . Here is one intermediate integral of order 2: 

st+l = 0. (13) 

The system f|T2l) + f|T3|) is involutive and so allows reduction to a rank 2 
distribution. Indeed this system can be re-written in the form similar 
to ([8]) — using Goursat representation ([9]) with ijj = — |A^/^ we getP 

^xx 3A ^ , Uxy A ^ , Uyy A ^ . (14) 

This system is compatible (compatibility writes formally as A^. = AA^ 
but this relation expresses the prolongation) and the characteristic field 



'Outside the submanifold in £^ given by Cartan equation 
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is ^ = Vx — Wy, which is also the Cauchy characteristic field for the 
Cartan distribution Cs on the equation [K4]. 

Quotient of this rank 3 distribution by ^ gives a rank 2 distribution 
on 5-manifold M^, which corresponds to the Monge equation 

w' = {z"Y/\ (15) 

This latter is equivalent to (I7j) and has Lie(G2) algebra of symmetries. 
Thus by Cartan's version of Lie-Backlund theorem |Ci| the contact 
symmetry algebra of the involutive system (fT2|) + (fT3|) is the same as 
for (ITSl) — it's one more representation of G2- 

Remark 2. These all are realizations of the non-compact (split) form 
of G2 as the maximal symmetric model. The compact form of G2 is 
realized by the automorphism group of the Calabi almost complex struc- 
ture {S^, J). By IK5I this is the maximal symmetric model that acts on 
non-linear overdetermined non-integrabl^ Cauchy- Riemann equations 

Higher Monge equations were studied in |AKj , and all maximal sym- 
metric models were identified as the following underdetermined ODEs: 
yM(a;) = (these again cannot be solved without quadratures 

according to Cartan's criterion). 

The symmetries here can be thought of as internal or external - they 
coincide by a version of Lie-Backlund theorem from [AKj . 



There are also PDE models for these according to [K4]. We will 
demonstrate some in examples, which also indicate a relation to the 
projective geometry of curves (the tool from [DZj ) — in our case (most 
symmetric models) these are the rational normal curves. 

Example 8. The following system is involutiv^ 

It has type SE's in notations of [K3J . Quotient by Cauchy characteristic 
= Vx — XDy yields a rank 2 distribution on a manifold of dimension 8. 
The weak growth vector of this distribution (we refer for the definition 
and properties to ( AK I|) is (2,1,2,3) and the corresponding Monge 
system is 

y"=\{z"')\ w' = \{z"'f. (IT) 

The corresponding graded nilpotent (Carnot-Tanaka) Lie algebra is 
free truncated of length 4 with 2-dimensional fundamental space g_i. 



^This means the corresponding Nijenhuis tensor Nj is non-degenerate. 

^It has Lie class w = 1, i.e. the solutions depend upon 1 function of 1 argument. 
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Similarly we reduce 4i?4 

Uxxxx 5"^ ' '^xxxy 4 A , Uxxyy 3 A , Uxyyy ' ^Wyy 

to a rank 2 distribution with growth (2, 1, 2, 3, 4) and the Monge equa- 
tion 

We can modify the symmetric models without destroying the sym- 
metry algebra. For the above SE^ we get its tangent 2D cone 

Uxxx = l"^^ + -^^/^j Uxxy = l"^^ + -^^Z^) 

'^xyy ~^ ^yyy A ~1~ /X. 

This type 2E^ system is compatible (the prolongation obeys A^.. = AAj^), 
and its general solution depends on = 2 functions of 1 argument. 

The characteristic is still — XDy (with multiplicity two) and the 
contact symmetry algebra is the same as the contact external algebra 
for (fT6l) and the same as the internal algebra for (fT71) : it has dimension 
12 and Levi decomposition g = SI2 ixt (radical r is 9-dimensional and 
it consists of 1-dimensional center and 8-dimensional nil-radical). 

Further modification gives the 3D tangent cone of the normal curve, 
which is a strictly parabolic 3rd order determined PDE with the same 
12-dimensional symmetry algebra q. Denoting by a,/3,7,5 the 3rd 
derivatives Uxxx,Uxxy,Uxyy,Uyyy we can write this system of equations 
parametrically 

a = iA"^ + AV + 3A^i/, 7 = iA^ + A/i-F z/, 
/3 = iA^ + AV + 2Az/, (5 = A + /i, 
or in the implicit form 

+ 807^(37 - 6^) + 27/3^27/3^6 - 8/36^ - 187^ + = 2187 /^Vs. 

Its solution space depends on w = 3 functions of 1 argument and is 
intrinsically related to the Monge system (fTTI) . 

Higher analogs of the above are valid on the basis of works [ AKl [K4 . 



4. Conclusion 

In this paper we discussed compatible overdetermined as well as 
underdetermined systems S of differential equations. We indicated that 
symmetries generically produce particular 'automodel' solutions to £. 
But if symmetries are few, they do not allow complete integration. 

On the other hand the abundance of symmetries is a sign of integra- 
bility. In Section |3] we briefly indicated via examples how large algebras 
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of symmetries make models unique. These symmetries often do not live 
on the equation-manifold, and a covering is required to see them. In 
this way Darboux integrability [Dl |F] was recast into the language of 
group quotients in [AFV] and new examples were integrated, see [AFj. 

Integration in closed form is too restrictive in the context of PDEs, 
and instead one considers reductions to ODEs. These latter are not 
arbitrary in the maximal symmetric they also possess sym- 

metries. This restricts the models (like those Monge equations from 
[AK] ) and gives a method to understand integrability. 
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